Mixed-valent regime of the two-channel Anderson impurity as a model for UBei3 
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We investigate the mixed-valent regime of a two-configuration Anderson impurity model for uranium 
ions, with separate quadrupolar and magnetic doublets. With a new Monte Carlo approach and 
the non-crossing approximation we find: (i) A non-Fermi-liquid fixed point with two-channel Kondo 
model critical behavior; (ii) Distinct energy scales for screening the low-lying and excited doublets; 
(iii) A semi-quantitative explanation of magnetic-susceptibility data for Ui-^Th^Beis assuming 60- 
70% quadrupolar doublet ground-state weight, supporting the quadrupolar-Kondo interpretation. 



PACS numbers: 75.20.Hr, 75.30. Mb, 71.27. +a 

Since the 1950's, Landau's Fermi-liquid theory has 
shaped our understanding of the metallic state. Based 
on the notion of a one-to-one mapping between the low- 
lying excitations of the interacting system and that of 
the noninteracting electron gas, the theory provides a re- 
markably robust scenario for the low-temperature prop- 
erties of interacting electron systems. It is against this 
outstanding success that a growing class of /-shell mate- 
rials — predominantly Ce- and U-based alloys — received 
considerable attention in recent years. Characterized by 
a logarithmically divergent linear coefficient of specific 
heat and anomalous temperature dependences of the re- 
sistivity and susceptibility, these materials appear to de- 
part from the conventional Fermi- liquid scenario Q , thus 
challenging our understanding of metallic behavior. 

In this paper, we present results on the two-channel 
Anderson impurity model in the mixed-valent regime, 
motivated by the unusual nonlinear susceptibility data of 
the non-Fermi-liquid (NFL) alloy system Ui-^Th^Beia. 
The restricted Hilbert space of our model, a ground 
quadrupolar (non-Kramers) doublet in the 5/ 2 config- 
uration and a ground magnetic (Kramers) doublet in the 
5/ 3 configuration, renders it intractable to study by con- 
ventional Monte Carlo methods. We have developed a 
new Monte Carlo method based upon the mapping onto 
a Coulomb gas. We use this method to calibrate non- 
crossing approximation (NCA) results, which can then 
be extended to more extreme parameter regimes. We 
find that the model displays NFL physics characteristic 
of the two-channel Kondo model, even at the extreme 
mixed-valent limit when the two charge configurations 
are degenerate. We also find that two energy scales ap- 
pear in the screening process away from the configura- 
tional degeneracy point. Using this model we are able to 
semi-quantitatively explain the linear and nonlinear sus- 
ceptibility data by assuming 60-70% ground-state weight 
for the quadrupolar doublet, with the Th doping induc- 
ing a higher 5/ 2 count, consistent with expectations from 



lattice constant data. We do not, however, explain the 
small energy scale for UBei3, and anticipate that dynam- 
ical inclusion of excited crystal field (CEF) levels in the 
5/ 2 configuration may remedy this problem. 

A proposed scenario for the NFL physics of UBei3 in- 
volves the screening of uranium quadrupole moments in 
the 5/ 2 configuration by conduction orbital motion [||. 
This quadrupolar Kondo effect can explain the enhanced 
specific heat and other data in this compound. In prin- 
ciple, though, the screened moment in the scenario of 
Ref. Q can either be magnetic or quadrupolar, depend- 
ing on which ionic charge configuration (U 3+ or U 4+ ) 
is lower in energy. Which case is realized in UBei3 is 
still unresolved. While the nonlinear magnetic suscepti- 
bility H is suggestive of a magnetic trivalent state, the 
rather weak temperature dependence of the linear suscep- 
tibility is more consistent with a nonmagnetic tetravalent 
state. Recently, following experiments on the nonlinear 
susceptibility of Uo.gTh .iBei3, Aliev et al. p suggested 
a different physical regime, in which strong quantum fluc- 
tuations drive the uranium ions to a mixed-valent state. 

Here we explore the possibility of a mixed-valent state 
in the framework of the two-channel Anderson model [0 , 
which consists of T% conduction electrons carrying both 
spin (a —] . I) and quadrupolar (a = ±) quantum labels, 
hybridizing via a matrix element V with a local uranium 
ion. The latter is modeled by a T 3 quadrupolar doublet 
in the 5/ 2 configuration and a Tg magnetic doublet in the 
5/ 3 configuration, separated in energy by tj = E(5f 2 ) — 
E(5f 3 ). The corresponding Hamiltonian reads 

W = E ^ c La^ a a + e f |5/ 2 , a) <5/ 2 , a\ (1) 

kaa a 

+ ^E{ c L.|5/ 2 ,-a)<5/V|+h.c.}, 

kaa 

where el creates a T$ conduction electron with spin a 
and quadrupolar moment a. This model may be applied 
to Ui_ x Th x Bei3 for temperatures above about half the 
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Kondo scale, for which coherent lattice effects are small. 

In the integer- valence limit, T = npV 2 <C |e/|, (p be- 
ing the conduction-electron density of states at the Fermi 
level) Eq. (|]) reduces to the two-channel Kondo Hamilto- 
nian Q, which is exactly solvable by a number of meth- 
ods: Bethe ansatz |Q, conformal field theory ||], and 
Abelian bosonization |9|. However, none of these solu- 
tions extend to the Hamiltonian of Eq. ([!]) , where charge 
fluctuations are present. This model is also intractable 
to conventional determinantal Quantum Monte Carlo. 

To overcome the difficulties associated with the Hamil- 
tonian of Eq. (Q) we devised a new approach, based 
on (i) mapping the corresponding partition function 
onto a classical one-dimensional Coulomb gas, and (ii) 
sampling the latter gas using Monte Carlo techniques. 
This approach allows an accurate calculation of the low- 
temperature thermodynamics of the model deep into the 
mixed-valent regime, all the way from weak to strong 
coupling. It also enables the computation of the nonlin- 
ear susceptibility, which, due to the high-order correla- 
tion function involved, is typically inaccessible to deter- 
minantal Quantum Monte Carlo. Below we outline our 
approach, starting with the mapping onto a Coulomb gas. 

The formal connection between the Kondo effect and 
the statistical-mechanical problem of a one-dimensional 
Coulomb gas was first recognized for the Kondo Hamil- 
tonian |[o|, and later extended to the one-channel An- 
derson model ■ Here we generalize the mapping to 
the Hamiltonian of Eq. (jl|) . The basic idea is to expand 
the partition function in powers of V, expressing it as a 
sum over all possible histories of the impurity. A history 
is a sequence of hopping events, or kinks, which overall 
preserve the impurity state and the occupation of each 
conduction-electron branch. Each history is represented 
by a sequence of impurity states, {70, . . . ,7„}, and a se- 
quence of imaginary times, {n, . . . , r„}, corresponding to 
the instances at which hopping events take place. Trac- 
ing over the conduction-electron degrees of freedom gen- 
erates a long-range "Coulomb" interaction between the 
kinks within a given history. These couple through the 
four-component charges e) = (8N^ + ,SNi + ,5N^,SNi^), 
where 5N aa is the change in occupancy of the er, a 
conduction-electron branch due to the i-th hopping event 
(i.e., the transition from state 7i-i to state 7^). Denoting 
the impurity- free partition function by Zq, one has 
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where E-y is the bare energy of the 7 impurity state (zero 
for 5/ 3 and e/ for 5/ 2 ); F(t) is the interaction strength 
between kinks with time separation r; and tq and r, i+ i 
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FIG. 1. The impurity susceptibility, x{T), for e/ = and 
different F/D. x(T) is defined as the response to a field that 
couples linearly to one of the impurity moments (for tf = 0, 
the magnetic and quadrupolar susceptibilities are identical). 
Error bars are smaller than the symbols used. For T < Tk, 
Tk the Kondo scale, x{T) has the logarithmic temperature 
dependence x{T) ~ (cl/Tk) log(Tjf /T), with a extracted from 
a logarithmic fit (dashed lines). In going from large to small 
F/D, a takes the values 0.019, 0.024, and 0.029, while T K /D 
equals 0.09, 0.032, and 0.00662, respectively. The solid lines in 
the middle and lower graphs are the results of the NCA with 
F/D = vr/23 and F/D — tt/144, respectively (see text). There 
is excellent agreement between the NCA and the Monte Carlo 
down to T /Tk ~ 0.7, at which point the NCA curves cross 
over to a logarithmic temperature dependence with a = 0.037 
and 0.042 for F/D = tt/23 and F/D = tt/144, respectively. 



are equal to zero and /?, respectively. At zero tempera- 
ture, F(t) reduces to ln(l+|r|/r c ), where D = l/r c plays 
the role of a bandwidth. At T > 0, the logarithm is re- 
placed by a more complicated expression. The neutrality 
condition fe^o selects only those histories that overall 
preserve the occupation of each conduction branch. 

The impurity contribution to thermodynamic quanti- 
ties can be computed directly from the Coulomb-gas rep- 
resentation of Eq. (||) using Monte Carlo. This formula- 
tion has the crucial advantage of being free of any sign 
problem, as the different terms in Eq. (||) are all positive. 
Below we present our results for the mixed-valent regime. 

We begin our discussion with the limit of strong va- 
lence fluctuations, |e/| <§C T, represented by e/ = 0. In 
Fig. 1 we have plotted the impurity susceptibility, x{T), 
in response to a field that couples linearly to one of the 
impurity moments — either the magnetic moment in the 
case of a magnetic field, or the quadrupolar moment in 
the case of an electric or strain field (the two suscepti- 
bilities are identical for e/ = 0). In a Fermi liquid, x(^) 
saturates at a constant as T — ► 0. Here it diverges loga- 
rithmically with decreasing temperature, consistent with 
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FIG. 2. Tk versus F, for e/ = 0. Open circles are results of 
the Monte Carlo; full line was obtained using the NCA. Error 
bars for the Monte Carlo data are smaller than the symbols 
used. Inset: The screening temperatures Tk /D (open circles) 
and T ex /D (filled squares) as a function of |e/|/D, for fixed 
T/D = 7r/16 (results obtained using Monte Carlo). 

the NFL behavior of the two-channel Kondo model [Jfj . 
The logarithmic temperature dependence extends from 
weak coupling (r/D as 0.03 in Fig. 1) to strong coupling 
(r/D ps 0.8). Thus, the same critical behavior of the 
integer- valcnt limit persists into the mixed- valent regime. 

The crossover to logarithmic temperature dependence 
in Fig. 1 is associated with a characteristic energy scale 
or Kondo temperature, Tk, defined as the tempera- 
ture at which the effective moment per unit occupancy, 
fi ef f(T) = Txi(T)/m(T), is 60% screened. Here m(T) 
and Xi(T) are the occupancy and susceptibility of the 
corresponding 5/' doublet [i = 2,3; we use ^b9i = 1, 
except in Fig. 3). Note that the /-electron count is given 
by 3 — ri2(T), and that rij(T) is fixed at one half for 
€f = 0. At high temperature, /x e //(T) reduces to the 
free-moment value of one quarter. Hence Tk is defined by 
/Lt e //(TV) = 0.1. For the cases shown in Fig. 1, this defi- 
nition gives good agreement (within 25%) with a slope of 
1/(20Tk) for the logarithmic component of Xi(T)/nj(T), 
as is characteristic of the two-channel Kondo effect . 

Figure 2 depicts the Kondo temperature Tk as a func- 
tion of r, for ej — 0. Rather than varying exponentially 
with 1/r, over the range 0.01 < T/D < 0.2 we find the 
power-law dependence Tk ~ T x , with x m 0.9. This 
behavior is very close to the linear dependence expected 
from simple rcnormalization-group (and analytic NCA) 
arguments. However, we find Tk/T m 0.2, which is not 
expected from the renormalization-group and NCA argu- 
ments. For a representative value T = 0.3eV, this ratio 
gives a Kondo temperature of Tk ~ 600K, i.e., sixty- 
fold larger than the 10K seen in UBei3. Hence, contrary 
to previous claims for an uranium ion with full spheri- 
cal symmetry fl4|| , our model does not support a small 
energy scale in the limit of strong valence fluctuations. 



Figure 2 also presents a comparison between the Monte 
Carlo approach and the NCA j|. Within the NCA, the 
low-energy physics is exclusively determined by the ra- 
tio of the spin degeneracy to the number of independent 
conduction-electron channels ||. Thus, irrespective of 
the inter-configurational energy e/, the NCA yields the 
critical behavior of the two-channel Kondo model for the 
Hamiltonian of Eq. (|l]), in agreement with the Monte 
Carlo result. In Fig. 2, the NCA curve for Tk is slightly 
shifted (on a log-log scale) with respect to the Monte 
Carlo one. This shift may be accounted for by rescaling 
T NCA relative to T MC , which we attribute in part to the 
different high-energy cut-off schemes used in the two for- 
mulations. By matching the NCA and the Monte Carlo 
occupation numbers away from ef = 0, we extracted the 
rescaling factor T MC /T NCA rj 1.44. Using this rescal- 
ing, one obtains good quantitative agreement between 
the Monte Carlo and the NCA throughout the mixed- 
valent regime, as is exemplified in Fig. 1 for the scaled 
susceptibility, Trtx(T). 

In going to e/ ^ 0, the spin and the quadrupolar mo- 
ments are no longer screened at the same temperature. 
Specifically, the moment associated with the excited dou- 
blet is quenched first at a characteristic temperature T ex , 
followed by the moment associated with the low-lying 
doublet which is screened at Tk < T ex . Physically, each 
screening temperature represents a different crossover. 
Screening of the excited doublet sets in as the occupancy 
of that doublet crosses over from a high-temperature free- 
ion (Boltzmann) form to one that is governed by valence 
fluctuations. On the other hand, Tk marks the crossover 
to strong coupling and the onset of NFL behavior. 

In the inset of Fig. 2 we have plotted T ex and Tk as 
a function of |e/|, for T/D = it/16. The ground-state 
occupancy of the 5/ 2 doublet ranges in this plot from 
n 2 = 0.29 for e f = 0.2 to n 2 = 0.71 for e f = -0.2. Note 
that T ex and Tk interchange meanings, depending on the 
sign of e/: for ef > 0, T ex and Tk are the screening tem- 
peratures for the quadrupolar and the magnetic moment, 
respectively; for e/ < 0, the roles are reversed. 

As seen in Fig. 2, T ex grows essentially linearly with 
\ef | . This stems from the fact that the free-ion occupancy 
of the excited doublet decays exponentially with |e/|/T. 
Consequently, the crossover to an occupancy driven by 
charge fluctuations is exponentially sensitive to \ef\/T ex . 
Contrary to T ex , Tk decreases with increasing |e/|, be- 
coming exponentially small for |e/| T. Extracting the 
ratio Tfc/r for |e/| = 0.2 and using the value T = 0.3eV 
one obtains Tk ~ 150K, which is still an order of magni- 
tude too large to explain the 10K observed in UBei3. 

Finally, we have computed the nonlinear susceptibil- 
ity! X^ 3 \ defined from expansion of the magnetization 
in the direction of the applied field: M(H) = x^H + 
■^X^H 3 + • • ■■ For UBei3, in addition to linear splitting 
of the 5/ 3 magnetic doublet, a finite magnetic field in- 
duces a van Vleck splitting and an overall shift of the 5/ 2 
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FIG. 3. x (3) (T) along the (100) direction, for F/D — it/ 16 
and tf/D = -0.2 (o), -0.08 (V),0 (A), and 0.2 (□). The 
corresponding ground-state 5/ 2 weights equal 0.71,0.59,0.5, 
and 0.29, respectively, while T K / D = 0.032,0.02, and 0.009 
for \ef\/D — 0,0.08, and 0.2. Inset: Linear susceptibility 
for the same set of parameters. For Tk ~ 10K, as appro- 
priate for UBei 3 , X (3) ( T = Tk) ~ -0.12emu/moleT 3 for 
ej/D = —0.08, in close agreement with experiment (cor- 
responding value for UBei3 is indicated by arrow). Similarly, 
X (1) (0) w 0.014emu/mole for the e f /D = -0.2 case, com- 
pared with an experimental value of 0.012 — 0.016emu/mole. 



quadrupolar doublet, due to matrix elements between the 
T3 ground doublet and excited CEF multiplets. We have 
included this effect by computing the T3 energy shifts 
perturbatively, up to fourth order in H. The first ex- 
cited CEF levels at Acef ~ 15mev |l5| were assumed 
to be degenerate T4 and triplets, chosen as to allow 
the van Vleck contribution to the linear susceptibility 
to quantitatively match the UBei3 data, as discussed in 
Ref. M. This also fixed the position of the Ti level. 
To make contact with experiment, we fixed the ratio 
Tk/Acef ~ 0.055, which is necessary to account for the 
overestimated Kondo temperatures in our calculations. 

The resulting x (3) {T) curves are shown in Fig. 3. From 
these calculations it is apparent that only for ef/D « 
—0.08, which corresponds to a ground-state f 2 occupancy 
of about 60%, can one obtain a curve quantitatively 
compatible with the pure UBei3 data. The curve for 
ef/D = —0.2 (about 70% ground-state / 2 weight) resem- 
bles that for Uo.gTho.iBei3. In contrast, the Cf/D = 
and 0.2 curves have too strong a temperature dependence 
to correspond with experiment, both for the nonlinear 
and linear susceptibility (inset to Fig. 3). In every case, 
the characteristic energy scale is too large for all of these 
mixed-valent runs to account for the experimental data. 
We conclude that the scenario of Aliev et al. [|| , suggest- 
ing that UBei3 is more strongly mixed valent than antic- 
ipated previously, and that doping with Th increases the 
5/ 2 weight (driving the system closer to the quadrupolar 
Kondo limit), is qualitatively consistent with our results. 



The most severe omission from our model is a proper 
dynamical treatment of excited CEF triplets (beyond per- 
turbative shifting of the bare levels). Although these 
triplets possess negligible Boltzmann weights in the ionic 
limit, they can have significant zero-temperature quan- 
tum weights for nonzero hybridization. Due to their mag- 
netic character, their inclusion can potentially render the 
%( 3 ) curves more negative even in the 5/ 2 limit. They 
also introduce (even in the mixed-valent limit) a small 
energy scale associated with crossover physics between 
ground-state screening of the T 3 doublet, and collective 



screening of the entire Hund's rule multiplet |16|. 
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